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Abstract. In this paper we establish a comparison result through symmetrization for so- 
lutions to some boundary value problems involving the fractional Laplacian. This allows to 
get sharp estimates for the solutions, obtained by comparing them with solutions of suitable 
radial problems. Furthermore, we use such result to prove a priori estimates for solutions in 
terms of the data, providing several regularity results which extend the well known ones for 
the classical Laplacian. 



1. Introduction and main results 

The final goal of this paper is to obtain a comparison principle using symmetrization tech- 
niques in order to get sharp estimates for solutions to some elliptic boundary value problems 
involving the fractional Laplacian operator. Let u : M. N — > R be in the Schwartz space of 
rapidly decaying functions, here the fractional Laplacian (— A) a / 2 of u with a G (0,2), is 
defined in a standard sense, that is either by the Riesz potential 



(-Ar/ 2 u(*) := C N , a P. V. / U(X) "ff g, 



where P.V. denotes the principal value and Cm, a is a suitable normalization constant, or as a 
pseudo differential operator through the Fourier transform on the Schwartz class 

j-[(-Ar/\](o :=wm{i) 

where J~[g] denotes the Fourier transform of a function g. As for the equivalence between 
these two notions, as well as a detailed description and properties concerning more general 
integro-differential operators, we refer to the book of Landkof [29J and the paper [37j . 
It is well known (see for instance [T7]) that to any function u smooth enough on 1^ we 
can associate its a-harmonic extension, namely a function w defined on the upper half-space 
:= M N x (0, +oo) which is a solution to the local (degenerate or singular) elliptic problem 

- div {y l - a Vw) = in 

(1) 

w(x, 0) = u(x) in Mr. 
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Moreover, Caffarelli and Silvetre give in |17) an interpretation of the fractional Laplacian 
(—A)"/ 2 as a Dirichlet to Neumann map: 

(-A) tt/2 ^) = 4 Km y- a ^(x,y), (2) 
k a y ^o+ ay 

where k a is a suitable constant, whose exact value is 

2 1 ~°T(1 — -) 

«« = r(f) 2 • (3) 

In order to define the fractional Laplacian in bounded domains the above characterization 
has to be suitably adapted. This has been done in the papers [12] and [14] . where formula 
Q allows to define the fractional Laplacian (— A) a//2 over a proper function space on fi, as 
we shall see in Section 2. 

The fractional Laplacian appears in several contexts. For instance, it arises in the study of 
various physical phenomena, where long-range or anomalous diffusions occur. Just to give few 
examples, this kind of operator can be found in combustion theory (see [19J), in dislocations 
processes of mechanical systems (see [25]) or in crystals (see [22]). Moreover, as it is well 
known in the theory of probability, the fractional Laplacian is the infinitesimal generator of a 
Levy process (see for instance [36]). Due to all of that, lots of authors devoted their interest 
to the subject. We just mention [37], [20], [IS], [H] dedicated to the obstacle problem and the 
free boundaries for the fractional Laplacian, the papers [H], |13j regarding some aspects of 
nonlinear equations involving fractional powers of the Laplacian, the convex-concave problem 
for the fractional Laplacian described in [T2], the work |27| in which a critical exponent 
problem for the half-Laplacian in an annulus is investigated, the study [30] of a nonlocal 
energy variational problem, and the papers [7J, [TUJ, [TT], [8]. Obviously this list is very far 
from being exhaustive. 

In order to describe our main result let us consider the nonlocal Dirichlet problem with 
homogeneous boundary condition 

'-A) a/2 u = f(x) infl 

(4) 

•i = on dQ, 

where O is an open bounded set of M and / is a smooth function on O. To this problem it 
is possible to associate the local one 

— div (y 1_a Vw) =0 in Cq 

w = on OlCq 



(5) 



lim y l a -^— = f (x) in Q, 

Ka y ^o+ dy 



where Cn := x (0, +00) is the cylinder of basis and OlCq := dQ x [0, +00) is its lateral 
boundary. We recall that (see for instance [12]) given a solution w in the weak sense to 
problem ([5]), then its trace on Q, trn(w) = w(-,0) =: u is a solution to problem Q (see also 
Section 3 for precise definitions). 
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Following [35], the idea is to get sharp estimates for the solution u to Q by comparing it 
with a solution eft to the radial problem 

{-A) a/2 <P = f#(x) inQ# 

(6) 

cf) = on dfi,^, 

where is the ball centered at 0, having the same measure as O and is the Schwarz 
rearrangement of /. Since <fi is the trace on of a solution v to the problem 

■ div (y 1_a Vu) =0 in 

= on d L C* ^ 

1 dv 

lim y 1_a 7^— = /* (aj) in 

Ka <,^o+ dy 

where := Q# x (0, +oo), <9lCq# := dVL# x [0, +oo), it makes sense to look for a comparison 
between concentrations of the functions w and v through their Schwarz rearrangements (see 
Section 2 for definitions). More precisely, we prove that 



ps ps 

/ w* (a, y) do < / v*(a,y)da Vs € [0, \Q\ 
Jo Jo 



(8) 



where w*(-,y), v*(-,y) are the one dimensional rearrangements of w,v respectively, for any 
fixed y 6 [0, +oo) . The achievement of such result looks reasonable because of the nature 
of problem for which a symmetrization with respect to x keeping the y variable fixed 
(i.e. Steiner symmetrization with respect to the line x = 0) is available. The key role in this 
framework is played by a second order derivation formula for functions defined by integrals, 
obtained in [T] for the smooth case and in |21j for less regular functions. 
We point out that through inequality ^ we easily get a comparison result between the traces 
on f2 x {0} of w and v, namely an integral comparison between u and 4>: 

Theorem 1.1. Let u and <fi be the weak solutions to problems Ml) and (|6|), respectively, and 

2N 

f G Ljv+«(J1), with a G (0,2). Then we have: 



/ u*(o)do< / (f>*(cr)da Vs G [0, |0| 
Jo Jo 



We emphasize that since for a = 2 the fractional Laplacian coincides with the classical 
Laplacian, for which comparison and regularity results via symmetrization methods are well 
known (see e.g. [38], [32], [2]), in the following we consider only the case < a < 2. 

In this setting the comparison result proved in Theorem 1.1 allows us to prove a priori 
estimates for solutions of problem Q in terms of the data /, providing several regularity 
results which extend the well known ones for the classical Laplacian (see Section 5). 
We also want to remark that the application of symmetrization techniques to general Levy 
processes is not new, as it is shown for example in [5] and [1] where several isoperimetric-type 
issues are investigated. Moreover, our approach is completely "PDE oriented" and it is not 
based on a probabilistic setting. 
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This paper is organized as follows. In Section 2 we give some basic definitions and proper- 
ties concerning the functional setting we are going to work with. In particular, in Subsection 
2.1 we recall some fundamental definitions concerning the fractional Laplacian in bounded 
domains and properties of weak solutions to the related Dirichlet problem. Subsections 2.2 
and 2.3 introduce the notion of Schwarz rearrangement and Lorentz space, together with some 
related properties. In Subsection 2.4 we define the Green function for the fractional Lapla- 
cian, an essential tool to write down an explicit integral representation of the solution <f> to the 
radial problem Q. In Section 3 we prove the comparison results stated above. Furthermore 
in Section 4, we exhibit some regularity results of the solution u in terms of the source data /. 
Finally in Section 5 we give some comments about the best constant in the L°° estimate. In 
particular we explicitly compute it for u on a ball, considering only the case a = 1 and N = 3. 



2. Preliminaries 

2.1. Function Spaces and Definitions. As we pointed out in the introduction, formula 
(2) given in [17] connects the nonlocal character of (— A) a ^ 2 to local problems of the form 
(1). This interpretation can be extended to the case of bounded domains. To this aim, it is 
convenient to introduce here a suitable functional setting and basic definitions. For all the 
details and proofs of the following definitions and properties, we refer to the the papers fT4"] , 
[12], [33] and gg. 

If is a bounded domain in M. N , the half-cylinder with base and its lateral boundary will 
be respectively denoted by 

C n := n x (0, +oo) and d L C n := dQ x [0, +oo). 

We introduce then the weighted energy space 



X£(C n ) := G H\C n ), w = on 8 L C n : jf 



y a \Vw(x,y)\ dxdy < oo 



equipped with the norm 



\w\\xg ■= ^jT y 1 a \Vw(x,y)\ 2 dxdy 



1/2 

y 1 ~ a \Vw(x, y)\ 2 dxdy 
IC n ) 

Thus we define the trace space by 

V a (fi) = {u = tr n w := u;(-, 0) : w G X£(C n )} , (9) 
where tr^ is the trace operator on the space w G Xq[Cq). 

The fractional Laplacian in a bounded domain 0, is well defined for function in V a (0). Indeed 
it is well known (see e.g. |14j . |12| ) that for any function u G V a (£l) there exists a unique 
minimizer uu to the problem 

infjy y l - a \Vw{x,y)\ 2 dxdy : w G X£(C n ), w{-, 0) = u in ftj . 
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By standard elliptic theory such minimizer w is smooth for y > and satisfies 

— div (y 1_a Vu;) =0 in Cq 

w = on d L C n (10) 

w(-,0) = u in ft. 

This yields to consider an extension operator in the following sense: 



Definition 2.1. Given a function u G V a (Q), the solution w to problem (10) will be said the 
a-harmonic extension of u on the cylinder Cn and will be denoted by Ext Q u. 

Then the fractional Laplacian operator can be defined through the Dirichlet to Neumann 
map as follows (see e.g. [2], [12"]): 

Definition 2.2. For any u G V Q (fi) we define the fractional Laplacian (— A) Q / 2 acting on u 
as the following limit (in the distributional sense) 

(-A) Q / 2 n(x) := -— ]imy- a ^(x,y), 
K a j/->0 Oy 

where w =Ext a (u) and K a is given by 

Let {ifk} be an orthonormal basis of L? (O) made by eigenfunctions of —A in $7 with zero 
Dirichlet boundary conditions and the corresponding Dirichlet eigenvalues. It is classical 
that the powers of a positive operator in a bounded domain, evaluated on certain function 
u, are defined through the spectral decomposition of u using the powers of the eigenvalues of 
the original operator. So in the case of the fractional Laplacian (—A)"/ 2 , if 



u = y^akPk 

k=l 

we must have 

(-A)"/ 2 n = 5> fc AfW (11) 

k=l 



This definition is coherent with Definition 2.2 since it is possible to give the following char- 
acterization of the trace space V a (C): 

Proposition 2.1. The space V Q (f2) defined in |$|) coincides with the space 

( oo oo ^| 

H := < u G L 2 ($7) | u = ^^ak^fk satisfying ^^flfcA^ 2 < oo > . (12) 



,2 W 2 

k=l k=l 



Moreover if u G V a (0) admits the decomposition u = Ylk^=i a k i Pk; then its a- harmonic has 
the following explicit representation 



Ext Q u (x, y) = y^afc^fc (x) p{>^ 2 y) (13) 

k=l 
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where p solves the problem 



P"(s) + 



a 



s 

,1 — OL J 



p'{s)=p{s) s>0 



lim y a p' (s) 



{ P(0) 



1. 



Therefore, using (13) and Definition 2.2, equality (11) easily follows. 
According to [12], we have the following definition of weak solution to problems of the type 
©: 

Definition 2.3. Let f G where a G (0,2). We say that w G X^(Cq) is the weak 

solution to problem ^ if for any test function tp G Xq(Cq) the following identity holds: 



,1-a 



Vw(x,y) ■ V(f(x,y) dxdy = K a / f(x) tp(x, 0)dx. 



(14) 



We note that for any test function ip G Xq(Cq), by the Sobolev trace inequality (see [T2] ) 

2N 

it follows that the trace <£>(-, 0) on Q, x {0} belongs to L N - a (Q), hence the integral at the 
right-hand side of (14) makes sense. Besides, the classical Lax Milgram theorem ensures that 
a unique weak solution w G Xq(Cq) to problem ^ exists. 

Then the definition of weak solution to problem Q is strictly related to the solution of ^ 
in the following sense: 

2N 

Definition 2.4. Let f G L N + a (0), where a G (0, 2). We say that u G H is the weak solution 
to Q if u = trQW, and w is the weak solution to problem 

We observe that if u is the weak solution to (|4j) , than its ct harmonic extension Ext a w is 
smooth for y > and decays to zero as y — > oo (see [T2]). 

Finally we point out that the space H defined in (12) is an interpolation space and it is 
possible to prove that (see [H] for the case a = 1 and [31], [28] for the general case) 

' H a / 2 {fl) if a G (0,1) 
// I H^\n) if a = 1 

H% /2 (n) if a G (1,2) 

where H a / 2 (£l) is the usual fractional Sobolev space, H* /2 (VL) is the closure of (Q) with 
respect to the norm ||-||#<*/2m) and 

Hl' 2 m := ju G ffVa ( n) : j[ ^dx < oo} , 
with d(x) := dist(x, dQ). 
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2.2. Basic facts about rearrangements. Let Q be a bounded open subset of M. N and u be 
a real measurable function on We will denote by |-| the iV- dimensional Lebesgue measure. 
We define the distribution function \i u of u as 

l_i u (t) = \{x £ n : \u (a?) | > t}\ ,t>0, 

and the decreasing rearrangement of u as 

u* (s) = sup{t > : fi u (t) > s} , s G (0, 

Furthermore, if cjat is the measure of the unit ball in l w and is the ball of M. N centered 
at the origin having the same Lebesgue measure as f2, the function 

u *(x) = v*{uj n \x\ n ) ,x£n# 

is called decreasing spherical rearrangement of u. For an exhaustive treatment of rearrange- 
ments we refer to [3], [26] and to the appendix of [39]. Here we just recall the well known 
Hardy-Littlewood inequality (see [21]) 

rM r 
u (x) v (x) I dx < / u* (s) v* (s) ds = / u^(x)v^(x)dx (15) 
Jo Jn# 

where u, v are measurable functions on 0. 

We point out that as we will deal with two variable functions of the type 

u:(x,y)£C n ^u{x,y)eR (16) 

defined on the cylinder Cn := x (0, +oo), measurable with respect to x, we can define the 
Steiner symmetrization of Cq with respect to the variable x, namely the set Cq := fi* x (0, +oo) . 
In addition, we will denote by fi u (t, y) and u* (s, y) the distribution function and the decreas- 
ing rearrangements of (16), with respect to x for y fixed, and we define the function 

u* (x,y) = u*(u N \x\ N ,y) 

which is the Steiner symmetrization of u, with respect to the line x = 0. Obviously is a 
spherically symmetric and decreasing function with respect to x for any fixed y. 

Now we recall two derivations formulas, that will turn out very useful in the proof of the 
main result. The following proposition can be found in [32], and it is a generalization of a 
well-known result by Bandle (see [3]). 

Proposition 2.2. Suppose that w G H (0, T; L 2 (Q)) for some T > 0. Then 

w* G fl' 1 (0,T;i; 2 (0,|n|)) 

and if \ {w(x,t) = w*(s,t)} | = for a.e. (s,t) G (0, x (0,T), the following derivation 
formula occurs 

dw f s dw* 

— (x,y)dx= -——(s,y)ds. (17) 

w{x,y)>w*(s,y) °V JO °V 

Moreover, what follows is a second order derivation formula due to Mercaldo and Ferone 
(see |21j), which is a suitable generalization of that contained in [1], where only analytic 
functions are considered. 
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Proposition 2.3. Let w £ W 2 '°° (Co,)- Then for almost every y G (0, +00) the following 
derivation formula holds: 

'" w (x,yf 2 



I 

Jw(x 



d 2 



W 



9 2 



(x,y)>w*(s,y) 



dy 2 



(x,y) dx = J ( a 



/ w* (a, y)da- 

JO Ju 



9~y 



-dU 



N-l 



w(x,y)=w* (s,y) 
\ 2 



+ 



dy 



w(x,y)=w*(s,y) 



\V x w\ 



\V x w\ 



w(x,y)=w* (s,y) 



\V X W 



dU"- 1 (x) 



2.3. Lorentz spaces. As we will deal with some sharp regularity results of the solution u to 
Q in terms of the data /, we introduce here basic notions regarding the functional spaces 
where / will be supposed to belong to. 

Let Q be a bounded open set of M. N . We say that a measurable function u : Q — > K belongs 
to the Lorentz L p,q (O) for < p, q < +00 if the quantity 



\ u \\LP'i(n) 



+00 



tPU*(t) 



sup tPU*(t) 
K 0<t<+oo 



< q < 00 
q = 00 



(18) 



is finite. We remark that for p > 1, and q > 1, the quantity in ( 18 ) can be equivalently defined 
replacing u* (t) with 

u **(t) = - [ t U *(s)ds. 



We stress that the L p,q — norm, for every 1 < p,q < +00, is rearrangement invariant, that is 



u 



\LP-i(Q) — \\ u ll£P.«(fi#V 

Besides, we emphasize that L p,q (CI) = L P (Q), L p,oc (CI) = A4 P (the Marcinkiewicz space) for 
any 1 < p < 00 and, for l<q<p<r<oo the following inclusion occurs: 

L°°(ft) C L r (Cl) C LP' 1 (CI) C L p ' q (Cl) C L p ' p (Cl) = L p (Cl) C L p ' r (Cl) C L p '°°(Cl) C L q (Cl). 

Now we recall a convolution inequality in Lorentz spaces due to O'Neil ([33]), which will be 
an essential tool to obtain some a priori estimates for solutions to problems of the type Q 
in terms of data belonging to Lorentz spaces (see Theorem 4.3): 

Theorem 2.1. Suppose that f G L Pl ' 91 (R N ) , g G L P2 '" 2 (R N ) where 



#1 



Then f*g£ L r ' s (R N ) where 
and t > 1 is any number such that 



1 1 

— + — > 1. 

Pi P2 



1 1 1 

- + --1 = -, 
Pi P2 S 

1 1 1 

~ + ~^ 7" 
qi <?2 t 



Moreover 



\\f*g\ 



< 3s 



LPi'n(R N ) \ \9\\lp2-12 
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2.4. Spectral decomposition of the solution. In this section we highlight some properties 
concerning the representation of the solution to the fractional Poisson equation by the Green 
function and the link with its spectral decomposition. According to what we have said in 
Subsection 2.1, it is always possible to get a spectral decomposition of the solution u to Q 
in terms of the Fourier coefficients of the source term /. Indeed, suppose that {tpk} is an 
orthonormal basis of I? (O) made by eigenfunctions of —A in ft with zero Dirichlet boundary 
conditions and the corresponding Dirichlet eigenvalues. Therefore, if u G H is the weak 
solution to problem Q, having the decomposition 



u 



(19) 



fe=i 



then the fractional Laplacian of u has the spectral decomposition (11). Thus if 



Cfc^fc, 



fc=i 



where Ck = (f, </?fc).L 2 (Q) are the Fourier coefficient of /, the Fourier coefficients of u are 

Ck 



ak 



X 



a/2 ' 



(20) 



Now, let us denote by Qr>{x,y) the Green function of a bounded domain D C 1^ for the 
fractional Laplacian (— A) a / 2 . Then we have (see |29| . [9], [6]) 

-(-A)«/ 2 g D (x,y) = 5(x-y) mV'(D). (21) 
Next, suppose that the function Qd has the following expansion, for any fixed y £ D: 

oo 

Go{x,y) = ^2ck(y)X k (x). 
fe=i 



Then equality (11) provides the following spectral decomposition for the fractional Laplacian 
of Q D : 



(-ArJ 2 g D (x,y) = J2K /2 Ck(y)Vk(x)- 

k=l 



(22) 



equation ( 22 ) links to 



If we multiply both sides of equation (21) by ip m and integrate over D with respect to x, 

-¥>m(y) 



oo „ 

y^A^ /2 c fc (y) / ip k (x)Lp m {x)da 



i.e. 



that is 



Cm{y) = - 
Qo{x,y) = - ^2 



a/2 



X 



<f k (x) ip k (y) 



k=l 



X 



a/2 



(23) 
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Hence from (19), (20) and (23) we easily infer that 



u 



\ a/2 Jn 

k=i H JD 



(y)dy = - GD(x,y)f(y) 

JD 



dy. 



(24) 



When D is a ball B(0,R), we shall frequently use the following explicit expression of the 
Green function (see [6], [29], [9]) 



SW^) = -2-*^^r\- 

where x, y £ B(0, R) and 



r (ir) r f*V 2 lx _ yr N 



S2' 



+1) 



N/2 



ds 



(25) 



(R 2 - \x\ 2 )(R 2 - \y 
\x — y\ 2 



We stress that (25) coincides with the Green function of classical Laplacian for a = 2. Clearly 
we have 

ab 



\yB(0,R) 

for -B(0, R) s. t. x ^ y, where 

' N\ p ^a\- 2 
^N/2 



x,y)\ < 



\x-y 



N- 



a := T 



, r (2 V ) r (2) R N-2 b 



oo af -l 



o (s + R 2 ) 



N/2 



ds. 



(26) 



(27) 



3. Comparison result 



The aim of this section is to obtain a comparison result between the solutions of problems 
([5]) and ([7]). The symmetrization method allows to obtain a priori estimates which are the 
main tools to obtain regularity results. 

Theorem 3.1. Let w and v be the the weak solutions to problems (|5| and (|7|), respectively, 

2N 

and f £ L N + a (Q), with a G (0,2). Then we have: 

w* (a, z) da < [ v*{a,z)da Vs € [0, (28) 



for any fixed z £ [0, +00) . 

Proof. We first observe that actually there is a clever way to rewrite equation in problem ^ , 
that is 

1 — a dw d 2 w 

A x w + — + — = 0. 

y dy dy 1 

As a matter of fact, if we follow [18J and make the change of variable 

a; 
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we find that problem ([5]) is equivalent to the Cauchy-Dirichlet problem 

in Cn 



' d 2 



w = 



on 8 l Cq, 



(29) 



^ -^(i,0) = K a a a - 1 /W into, 



where (3 := 2 (a — 1) /a. The aim is to compare problem (29) with the corresponding sym- 
metrized one: 

,d 2 v 



z^ + A x v = 



in C$ 



v = 

dv 
dz 



on d L C* 



(30) 



Now we recall that w is smooth for any z > 0, so if for a fixed z > we consider the test 
function 



sign (w(x, z)) 
w(x, z)\ — t 



h 







if \w(x, z)\ >t + h 

sign (w(x, z)) if t < \w(x, z)\ <t + h 
if \w(x,z)\ < t, 



we can multiply the first equation in (29 ) by ip y h (x) and integrate over Q. A simple integration 
by parts yields the identity 



IVt-ioI dx — z' 1 



t<\w\<t+h 



\w\>t+h 



d 2 w p l 
dz^ dx - Z h 



d 2 w ( \w\ — t 



t<\w\<t+h 



dz 2 \ h 



sign (w) \ dx = 



Letting h — > and using the isoperimetric inequality, by standard arguments (see e.g. |38j)we 
get 



J a 



w(x,z)>t 



d 2 w 
dz 2 



dx 



dt 



N 2 cj» (fi w (t)) 2 -* < 0. 



Now if we set 



U (s, z) = I w* (a, z) da, 



using the second order derivation formula of Proposition 2.3 we find that U verifies the 
following differential inequality 



8 2 U 



d 2 u 



dz 2 



(31) 
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for a.e. s £ (0, |f2|) and for any z G (0, +oo) , where p (s) = N 2 ouj^ s 2 ~n . Moreover, the first 
order derivation formula (17) implies 

dU d F *, \i 9 f / f dw , ^ J 

-j— = w (a, z) da = — w(x,z)dx = — (x,z)dx, 

OZ OZ Jo OZ Ju)(x,z)>w*{s,z) Jw(x,z)>w*(s,z) c)z 



hence making use of the Hardy-Littlewood inequality (15), we easily get 

^(s,0) = / ^ (x,0)dx = -a a ~ l K a [ f(x)dx 

oz Jw(x,0)>w*(s,0) ° z Ju(x)>u*(s) 



>-a a ~ l K a f f*(a)da, s € (0, \Q\ 
Jo 



So the function U satisfies the following boundary conditions 
U(0,z) = Vzg[0,+oo) 

dU 

— (\n\,z) = VzG[0,+oo) 



qtj r s 

{s^)>-a a ~ l K a f*(a)da, s € (0, . 
Jo 



dz 



Now if v is the solution of the symmetrized problem (30), being v radially decreasing with 
respect to x, we obtain 

a d 2 V d 2 V 

P«^rir = (32) 



dz 2 ds 2 

where 



V (s, z) = v* (a, z) da. 
Jo 

Concerning the boundary conditions, we remark that in this case one has 
^(s,0) = - a ^ Ka [ f#{x)dx 

az J v(\x\)>v* (s) 

r(s/uj N )y N 

= -NuiNa*- 1 ^ / /* (oj N r N ) r N ' l dr 



K a [ f* (a) da s e (0,|O| 
Jo 



therefore V satisfies the conditions 

V(0,z) = VzG[0,+oo) 

dV 

— (\Q\,z) = We[0,+oc) 



If we put 



OT t 

_ ( S) o) = -a- 1 ^ J Q S f* (a) da, s 6 (0, . 
Z (s, z) = U (s, z) — V (s, z) = / [w* (a, z) — v* (a, z)]da 
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by (31) and (32), one has 



for a.e. (s,z) G D := (0, |0|) x (0, +oo) and the following boundary conditions hold 

Z(0,z) = Vz€[0,+oo) 

dZ 



In particular 



dZ 



ds 
dz 
>,0) 



(|J2|,z) = Vzg[0,+oo) 
(s,0)>0 s€(0,|fi|). 
BZ 

(s,o)<o s€(o,|n|), 



(33) 



ft,—' ft,—-" ~™ (34) 
where f is the outward normal to the line segment (0, \ We observe that the operator L 
is elliptic in any point (s,z) £ D hence by Hopf's maximum principle (see [35]), Z attains its 
maximum on the boundary of D, and in the points where the maximum is attained we get 

dZ 



du 



> 0. 



Hence by ( 33 ) , ( 34 ) , this ensures that 

Z(s,z)<0 s€[0,|fi|] 

that is 

rs rs 

/ w* (a, z) da < / v* (a, z) da s € [0, \Q\ 
Jo Jo 

for any z £ [0, +oo) . 



□ 



Obviously, since 4> the trace on of the solution v of ([7]) and u the trace on O, of the 
solution it; of pL by Theorem 3.1 we get Theorem 1.1. 



4. Regularity results 



In this section we are interested in regularity results for solution u of problem Q. Using 
Theorems 1.1 and 3.1, we are able to prove some regularity results of the solution u in terms 



of the data /. In the following we will use the integral form (24) for the solution <f> to the 
symmetrized problem Q, namely 



<j>(x) 



G n #(x,y) f*(y) dy. 



(35) 



Theorem 4.1. Let u the solution to problem (J4J), where f £ La ^(O.) with < a < 2. Then 

ueL°°(n). 



We start by generalizing a well-known result for the classical Laplacian: 
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Proof. Let us consider the solution <f> to problem Q. Since (p is radially decreasing, using 
(35) and (26) we obtain that, for some constant C, 



L°°(f7#) 



0(0)=/ f*(y)\g n# (0,y)\dy<*b f l^dy 
Jn# Jn# \y\ 



fny) 
\y\ N - 

Rn rM 



ab / r a - 1 f*(u; N r N )dr = a b I s^f(s) ds = a b||/|| v , 



On the other hand, Theorem 1.1 gives 

IMU°°(n) < ||0||ioo ( n#) 

and the result follows. □ 

Remark 4.1. We stress that if f G -L p (£l), /or some p > iV/a, then according to Lorentz 
embedding (see Subsection 2.3) by Theorem 4.1 we get u G L°°(0). 



A consequence of the comparison result of Theorem 3.1 is the boundedness of the a- 

— JV i 

extension w of u in C n when / G L"^"' for < a < 2. To prove this result, we first 
compute the solution v to the radial problem (30) by using the separation of variable method. 
We look for a function v, radial with respect to x, such that 

v(x,z) = X(\x\)W(z). 

Putting v inside the first equation of p0| ), we find there must be a value A such that 

-A X X = A (36) 




W{z) 

that is the function X{x) = X{\x\) solves the classical eigenvalue problem for the Laplacian 

in 

(37) 

on dft^, 
while M^(z) verifies the problem 

zPW (z) - XW (z) = 

(38) 

lim W (z) = 0. 

z— >+oo 

Therefore (X,X) = (Xk,Xk), for some k, where {Xk} and {Xfc(|x|)} are the eigenvalues and 
the radial eigenfunctions of the Laplace operator in 0* with zero Dirichlet boundary values 
on 30*, namely 

Xk= {lt) fc=1 > 2 >--- ( 39 ) 

where 



Rn 
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is the radius of the ball k are the zeros of the Bessel function J(n-2)/2( z ) °f order 
(N -2)/2., and 



1 / 2 \ 1/2 

Xk(r) = Rn\JN(O k )\ \N^J 



_ N-2 
r 2 Jat-2 



6k 
Rn 



A; = 1,2, 



(40) 



where r := \x\. We recall that the system forms an orthonormal basis of the space 

L^. ad (Q#) made by all radial functions in L 2 . 

Then, recalling that the solution <f> of Q is radially decreasing, we can represent it by 



^a fc X fc (r), 



(41) 



fc=i 



where the a k are given by ( 20 ) , and c k are the Fourier coefficients of with respect to ( 40 ) , 
i.e. 

--Rn 



c fc = Noo N I r"' 1 f*(u N r")X k (r)dr. 







Now, to each eigenvalue X k we associate a solution W k to problem (38). The equation in 
( 38 ) is a modified Bessel equation (see (Mj , [30] ) , whose solutions are combinations of Bessel 
functions of the third kind. According to the asymptotic behavior at infinity of the Bessel 
functions ([30]), we have that 

W k (z) = C k H k (z), (42) 

where 



H k (z) :=^~zK^_ ( 

2-/3 \ 



2-/3 



\ k Z 2 



(5 = 2 (a — 1) /a, the C k are constants and K u (t) is a Bessel function of the third kind. We 
also notice that 



H' k (z) 



2-/3 \2 — p 



AfcZ 2 



(43) 



Finally, using the boundary condition of problem (30), we can write the following explicit 
expression of v (here r = \x\): 



v{r,z) = Y,Xk{r)W k { 
k=l 

with coefficients 



1 



1/2 oo „ 



0k 

R~n 



C k H' k (0) 



(2iVw 7 v) 1/2 a Q " 1 K, 
Rn 1^(^)1 



JV 

r 2 Jjv_2 
2 



^ ] r^NT N )dr. 



RnJ 



r) H k (z), (44) 



(45) 



Of course the trace v (r, 0) given in (j44J) coincides with the solution represented by (41). 
Indeed by the asymptotic behavior (see [30] ) 



(46) 



then by (43), (45) and (41) we find 
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(r,0) = f2 X k(.r)C k H k (0) = 

k=l 

= ( 2N ^ 1/2 f [ Rn t m (h. t ) r 

Rn kTi \ Jn (e k ) Jo 2 \ R n J 



(u> N t N )dt 



Now we are able to prove the following result. 

Theorem 4.2. Let w the solution to problem where f G L"^ ,1 (il) with < a < 2. Then 
w G L°°(C n ). 

Proof. Let v be the solution of ([7]) as in (44). By the asymptotic behavior of the Bessel 
functions K u at infinity (see [30]) we deduce that H k (z) — > as z — > oo, therefore u(r, z) — > as 
z — > oo. Besides, since v(x,0) = </>(x), by Theorem 4.1 we find that v G L°°(Cq#). Moreover, 
Theorem 13.11 assures that 

lk(->*)IU«(n) < ll u (-^)llx°°(n#) Vzg[0,+oo), (47) 

hence u> G L°°(Cq). □ 



We emphasize that the result of Theorem 4.2 is not new (see for instance [12], [13]), al- 
though our techniques make us able to achieve the sharper L°° estimate (47). 

Now we provide new regularity results when / belongs to Lorentz spaces L(p, r) for p < N/a, 
obtaining the generalization of the corresponding classical regularity result for the Laplacian. 

Theorem 4.3. Let u be the solution to problem where f G L P)T {£1) with 



2N N 

— < p < — 

Jy + a a 



and r > 1. Then u G L q ' r (Q) with 



q :-- 



Np 



N — ap 

Proof. Extending / to zero outside f2*, inserting inequality (26) into (35) we find 

\<t>(x)\ <ab(f**\x\ a - N ). 



Then applying Theorem 2.3 with the choices g 
q2 = oo, we have s = q = Np/ (N — ap), t = r and 



■r 



a-N 



Pi = P, P2 = N/(N - a), qi = r, 



Li> r (yi*) 



< ab 



/* * \x 



a-N 



L9. r (n#) 



< 3gab||/|| L p,r {n#) |||a: 
= 3gab||/|| LP , r(n #). 



\a-N\ 



\ L N/(N- a ), 



Finally by Theorem 1.1 we get 



M|£9,r(n) 



< 



and inequality (48) allows to conclude. 



(48) 



□ 
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In virtue of Theorem 4.1 if / € L N ^ a ^(fi) there is a constant C such that 



(49) 



Due to the form of the Green function in ( 25 ) , it seems quite difficult to face the problem of 
finding the best value for C in (49). Nevertheless, we remark that this becomes reasonably 
easy when one replaces the L N l a ^- norm of / at the right-hand side of (49) with an IP norm, 
for some p > N/a (which is possible, by Remark 4.1). In fact, since the solution (j> is radially 
decreasing, in order to get an L°° estimate of 4> it is enough to look for a sharp upper bound 

(50) 



of (j)(0). To this end, we first observe that (35) yields 



where 



0(0) 



r(f)_ r /n 



^((s/uN^r^ds, 



t 



a-N / ' 



S2 







+ 1 



N/2 



ds. 



(51) 



We remark that it is possible to write an explicit form of the integral at the right-hand side 
of (51). Indeed, we know that (e.g. see |23j) 



S2" 







+ l) N /'- 



ds = 2-^ 2 iq 



a 



N a a ^ 2 



where 2-^1 (• ; • ; • i z) denotes the Gauss hypergeometric function. Therefore by (51) 

n\*/ N 



iP((s/u N ) 1 / N )=B N , a s^ 
where 

&N,a '■ 

So if we set 



^(l^-^)" /2 ,n(f; 

2 1 -°r(f) 



N a a 



(>/iV _ M 2/N 



3a -2 
N 



l-(a/N) 
N 



a/2 



2F1 



\n\*/ N 



N a a 

rr l+ r^/N s 2/ N 



using Holder inequality in ( 50 ) we have 



< \B 



N,a\ 



1J> (SI) 



\Q\ 



s^p' [ip(s)] p ' ds 



l/p> 



We point out that the function 92 is bounded in [0, \ (see also the picture below), so the 
integral at the right-hand side of the last inequality converges if and only if p > N/a. 
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Case : N = 3, a = 1, Q = B 1 



10 



The best constant C(N,p,a,Q) in (49) is then 



C(N,p,a,n) := | B N . 



|fi| 



i/p> 



o 



Example 5.1. Let us calculate the best constant C in the case of the square root of the 
Laplacian \J— A ("i.e. £/ie case a = 1), when N = 3 and $7 = £>(0, 1). In this case, we have 
the following, explicit form of the Gauss hypergeometric function: 



3 13 
2' 2' 2' 



2-^1 1 t:; tt; o'- 2 



i 



Then 



47T 



1/3 



4 \ 2/3 

3 



and we aaue, fry a change of variable, 

(2tt) 1 /p' 



C(p) 



2vr 2 



i -, , \W 

t*- p ' (1 - t)^ dt 



o 



(2 7 r) 1 / p ' / p-3 3p-2 
B 



2tt 2 V 2 (P- !)' 2(p- 1) 
where B(-,-) is the Euler beta function. 



(p-i)/p 
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